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Abstract 

It is shown how simple assumptions lead to piecewise linear behavior, which is observed 
in certain phase transitions. 

In certain phase transitions, some property F(cp) of a system is a piecewise linear function of a 
parameter (p. Its slope s(cp) = dF((p)/d(p changes discontinuously at a critical value tp c of cp. An 
example is provided by the force F(tp) required to push a rod into a layer of glass beads of volume 
fraction if, .57 < ip < .63, at constant velocity (Schroter et al., [HI). When the penetration depth is 
60mm, F{ip) is nearly piecewise linear, with its slope changing nearly discontinuously at (p c = .60. 
( See Figure 1.) We shall show that simple assumptions lead to piecewise linear behavior of F(cp). 
Similar assumptions lead to the sharp transition from one power law to another power law, which 
occurs in turbulent boundary layer flow. 0. 
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Suppose the derivative of the slope s(ip) is a smooth function f(s, ip), which we expand in a 
Taylor series about some value s , through quadratic terms. Then 

= V 9 ) = f( S 0, V) + fa(s , ip){s - S ) + ^fss{s , (f){s - S ) 2 + O [(s - S f] . (1) 

We write the quadratic polynomial in s, on the right side of (QQh in terms of its roots a and (3. We 
assume them to be real, with a < (3: 

^ = 1 (s-a)((3-s) + O[(s-s ) 3 ]. (2) 

When we omit the error term in ©, and approximate a, (3 and 7 by constants, we can solve the 

resulting differential equation for s((p). Upon setting s(ip c ) = (a + (3)/2, we get 

a + p e (P-<*h(v-<Pc) 
S W) = 1 + e (/*-<*h(¥>-¥>c) • (3) 

As 7 — > 00, the slope s(ip) given by © becomes discontinuous at ip c , tending to a for ip < ip c and 
to (3 for (p > ip c . For 7 finite but large, s(ip) changes continously but rapidly from a to (3. 
F((f) is obtained by integrating the defining equation dF/dcp = s(ip), with s given by ©: 



F(<p) = F(ip c ) + f s{y'W 

J(p c 



= F{ip c ) + p(<p -<p c ) + - [log (1 + e 0>-"M*-^>) _ log 2 ] . (4 ) 
As 7 — > 00, F((p) given by © tends to the piecewise linear function 

lim F(cp) = F((p c ) + a((p-(p c ), ip < (p c , 

7^00 

= F(ip c ) + (3(ip - ip c ), <p>ipc- (5) 

This is the form of F(ip) shown in Figure 1 and in Figure 2 of HI. For 7 finite, © provides a 
smooth approximation to the piecewise linear function in ©. 



In [|2| we treated the mean velocity u(y) at distance y from a wall in a turbulent boundary layer 
flow. We defined the dimensionless shear s(y) = d\ogu/d\ogy, and proceeded as in (QQ) - © 
with (p — y and (p c = y c . Then we integrated the defining equation d log u/d logy = s(y), with 
s(y) given by ©, to obtain u(y). As 7 — > 00, the result © for s(y) tends to a for y < y c and 
to (3 for y > y c . The result for u(y) tends to Ay Q for y < y c and to B?/ 3 for y > y c , with A 
and B constants. Previously these two power laws had been shown to fit the experimental data by 
Barenblatt, Chorin, and Prostokoshin [f3l. 
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Figure 1. F(<p) versus cp given by © with (p c = .598, F(cp c ) = .5, a = 10.7 and (3 = 80. 
These values were chosen to make the graph resemble that in Figure 2 of flU, which shows the 
force F((p) in newtons versus the volume fraction (p for a penetration depth of 60mm. 
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